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ABSTRACT 

Context. The vortex coronagraph is an optical instrument that precisely removes on-axis starlight allowing for high contrast imaging at 
small angular separation from the star, thereby providing a crucial capability for direct detection and characterization of exoplanets and 
circumstellar disks. Telescopes with aperture obstructions, such as secondary mirrors and spider support structures, require advanced 
coronagraph designs to provide adequate starlight suppression. 

Aims. We introduce a phase-only Lyot-plane optic to the vortex coronagraph that offers improved contrast performance on telescopes 
with complicated apertures. Potential solutions for the European Extremely Large Telescope (E-ELT) are described and compared. 
Methods. Adding a Lyot-plane phase mask relocates residual starlight away from a region of the image plane thereby reducing stellar 
noise and improving sensitivity to off-axis companions. The phase mask is calculated using an iterative phase retrieval algorithm. 
Results. Numerically, we achieve a contrast on the order of 10"^ for a companion with angular displacement as small as 4 A/D with 
an E-ELT type aperture. Even in the presence of aberrations, improved performance is expected compared to either a conventional 
vortex coronagraph or optimized pupil plane phase element alone. 

Key words, instrumentation: high angular resolution - planets and satellites: detection 


1. Introduction 


The vortex coronagraph (VC) is an optical system for high- 
contrast imagingofastronomicalo^ect^t small angular sepa¬ 
rations ( [Mawet et al.|2005t|Foo et aL|2005| ). The VC suppresses 
the light from a star allowing direct detection of dim compan¬ 
ions, exoplanets, and circumstellar disks. Imaging objects with 
a VC that are otherwise buried in the noise associated with the 
bright host star has been demonstrated in laboratory (e.g. |Mawe^ 
et al.|2009 Delacroix et al.|2013) and on-sky observations (e.g. 


Swartzlander et al.|2008HMawet et al.|2010[|Serabyn et al.|2010 


Absil et al.|2013[|Defrere et al.|2014|) 

The VC was originally devised for telescopes with a circu¬ 
lar aperture, where light from an on-axis point source is com¬ 
pletely rejected while light from off-axis sources is preserved 
( [Mawet et~ar|2005[|Foo et aL]|2005| ). Unfortunately, most tele¬ 
scope pupils have an obstructing secondary mirror with spider 
support structures, for which the on-axis source is only partially 
suppressed (e.g. |Jenkins||2008| ). Moreover, very large apertures 
are often formed by a segmented mirror, which is generally not 
circular and may have discontinuities between segments. 

Several solutions have been proposed to improve contrast 
performance of a VC on telescopes with complicated apertures. 
The effect of the secondary mirror may be mitigated by use of a 
sub-aperture ( [Mawet et al.|2010[|Ruane & Swartzlander|2013|, 
tandem coronagraphic stages ( [Mawet et al.|2011[[Gahcher et al. 
Mawet et al.[[2013a| ), or a ring-shaped apodizer ( [Mawet 


2011 


et al. 


2013b[ ). Spiders and other aperture disconti nuities may be 


compensated for by binary amplitude apodizers ( [Carlotti et al 


F.R.S.-FNRS Research Associate 


2014| ), focal plane phase co rrections ([Ruane et al.|2015a|b| ), or j 


pair of deformable mirrors ( [Pueyo & Norman|2013[ ). 

Whereas a VC uses a focal-plane phase mask and down¬ 
stream aperture stop ("Lyot stop") to suppress starlight, an al¬ 
ternate class of pupil-only coronagraphs have achieved consid¬ 
erable success on telescopes with complicated apertures, using 
amplitude ([Kasdin et al. [[2003 [ [Carlotti et al. [201 1[) and phase 
( [Codona & Angel[2004[ [Ken worthy et al.|2007[[2010[ ) pupil ele¬ 
ments. 

Here, we present continuous pupil-plane phase elements that 
introduce a spatially variant phase shift in the plane of the Lyot 
stop (see Fig. [^, which acts to relocate residual starlight away 
from a dehned region of interest where dim sources may be de¬ 
tected. The Lyot-plane phase mask (LPM) is intrinsically loss¬ 
less, may be designed for a large variety of apertures, and is sim¬ 
ple to integrate with conventional coronagraph designs. We show 
that improved contrast performance may be achieved compared 
to pupil-only phase mask coronagraphs. 


2. The vortex coronagraph 

In this section, we demonstrate that the contrast performance 
of a conventional VC is limited on telescopes with complicated 
apertures. But hrst, we briefly review the case where the pupil 
is circular and has no obstructions. The layout of the VC is il¬ 
lustrated in Fig. A vortex phase mask (VPM) is placed in 
the focal plane (FPl) of a 4-/ optical system with transmission 
t = exp(//0), where 0 is the azimuthal angle in FPl and I is 
an integer known as the "topological charge." A nonzero, even 
value of I is required for ideal starlight suppression with a cir- 
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Fig. 1. Schematic of a vortex coronagraph. The first pupil plane (PPl) is 
typically an image of the telescope aperture. Lens LI forms focal plane 
FPl with field F(p, 0). A phase-only element known as a vortex phase 
mask (VPM) is located at FPl, with transmission t - exp(//0), where / is 
a nonzero, even integer. Lens L2 forms the output pupil plane (PP2) with 
field E(r, 0). The on-axis starlight {a - 0, red rays) is diffracted outside 
of the Lyot stop (LS). Lens L3 forms the subsequent focal plane (FP2) 
with the on-axis starlight removed, while light from off-axis sources 
{a ^ 0, blue rays) propagates through the LS. A Lyot-plane phase mask 
(LPM) may be introduced at PP2 to improve the contrast between off- 
axis and on-axis sources in the image at FP2. Lenses LI and L2 have 
focal length /. 


cular aperture. Various techniques are available to fabricate such 
achro matic phase masks ([Bomzon et al.||2001[ [Marrucci et al 


2006[ [Murakami et al.|2012| ). The scalar field immediately after 


the VPM owing to an on-axis point source may be written 


^(P, 0) 


ka^ Ji (kapjf) 
/ kapjf 


( 1 ) 


where p is the radial coordinate in FPl, a is the pupil radius, 
k = 2711 A, A is the wavelength, / is the focal length, and Jn is the 
^th order Bessel function of the first kind. The field at the output 
pupil (PP2) is given by the Fourier transform of Eq. 

oo 

kcL C* 

E (r, 0) = y J {kapjf) Ji {krpjf) dp, (2) 

0 


where (r, 0) are the polar coordinates in PP2. Remarkably, Eq.[^ 
evaluates to a discontinuous function, which is zero-valued 
within the geometric image of the pupil (r < a) for nonzero, 
even values of /; that is, a nodal area appears. All of the light 
from a distant on-axis point source appears outside of this re¬ 
gion; for example, E (r, 6) = (alr)^e^^^ for r > a and I = 2. A 
circular aperture stop known as the Lyot stop (LS), with radius 
Gl where Gl < g, is placed in PP2 to block all of the light from 
the on-axis source. Off-axis sources (i.e. a 0) do not form a 
nodal area and therefore transmit through the LS. For a AjD, 
the transmitted energy increases as . 

Although the VC theoretically provides ideal starlight sup¬ 
pression with an unobstructed circular pupil, telescope apertures 
are often more complicated. For this discussion, we consider the 
future European Extremely Large Telescope (E-ELT). Figure 
shows the entrance pupil of the E-ELT (see Fig.j^a)) and a pos¬ 
sible LS (see Fig. |^b)) along with the complex fields at PP2 
directly after the LS owing to an on-axis point source and a VC 
with I = 2 (VC2, see Fig. ^c)-(d)) and / = 4 (VC4, see Fig. 
[^e)-(f)). The fields here, and throughout this work, are calcu¬ 
lated using the Fast Fourier Transform (FFT) algorithm on large 
computational arrays (16384 x 16384 samples), with 1046 sam¬ 
ples across the pupil diameter D and 15.7 samples per AF# in 
the image plane, where F# = fjD. It can be seen in Figs.j^c) 
and [^e) that light from the on-axis point source leaks through 
the LS. We note that the LS is slightly smaller than the entrance 
pupil with a dilated central obscuration and spiders reducing the 
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Fig. 2. (a) Pupil of the E-ELT. (b) Lyot stop designed for a VC on E- 
ELT. (c) Amplitude and (d) phase of the field directly after the LS for 
an on-axis point source and / = 2. (e)-(f) Same as (c)-(d), but for / = 4. 


maximum system transmission (i.e. without a focal plane ele¬ 
ment) by 4%. 


3. Lyot-plane phase masks for E-ELT 

In the case of the E-ELT pupil, the VPM and LS (together, the 
VC) block 95.5% of the total power due to the star with 1 = 2 
and 95.3% with I = 4. The residual starlight reaches FP2 where 
it is spread over several units of TE#, along with its associated 
noise (see Fig.[^a)-(b)). Improved sensitivity in a given discov¬ 
ery region may be achieved by designing a dark hole in the point 
spread function (PSF) of an on-axis source (see Fig. [^c)-(h)). 
This approach has found success in the apodizing phase plate 
(APP ) coron agraph ( jCodona & Angel|[2004t [Ken worthy et al.j 
[2007 1 12010| ). In comparison to a focal plane coronagraph, the 
APP uses a pupil-plane phase mask to produce a dark hole in 
the spatially invariant PSF rather than reducing the amount of 
on-axis starlight reaching the image plane. Our approach com¬ 
bines focal plane and pupil plane coronagraphy to first suppress 
the star with a VC and then sculpt a dark hole in the residual 
starlight with an LPM, further improving the contrast locally. 

To reduce unwanted starlight, we present two LPM designs 
for VC2 (LPMl and LPM2, see Fig. |5a)-(b)) and VC4 (LPM3 
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Fig. 3. Monochromatic on-axis PSFs (in log irradiance) for (a)-(b) VCs, 
(c)-(f) VCs with LPMs, and (g)-(h) EPMs only. Examples phase masks 
that form semi-annular and full-annular dark regions are shown. Each 
is normalized by the peak value of the E-ELT PSE. 


and LPM4, see Fig. [2c)-(d)). LPMl and LPMS form a semi- 
annular dark hole in FP2 (see Fig. [3 c),(e)), while LPM2 and 
LPM4 clear the full annulus (see Fig. ^d),(f)). We also present 
entrance pupil masks (EPM) that form semi and full annular dark 
holes (respectively EPMl and EPM2), without the use of a VPM 
in EPl (see Eig. 3(g)-(h) and Eig. 4 (e)-(f)). The EPM is sim ilar 
in principle to an APP coronagraph ( [Codona & Angel|2004| ). 

The LPMs and EPMs are calculated using a point-by-point 
iterative phase retrieval method ( [Ruane et al.||2015a|b| ), where 
the star is assumed to be an infinitely distant point source and the 
field in PP2 is found computationally using paraxial Eourier op¬ 
tics methods. The phase in PP2 is optimized such that a dark hole 
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Fig. 4. Optimized pupil plane phase masks for (a)-(b) VC2, (c)-(d) VC4, 
and (e)-(f) no VPM. The corresponding PSFs are shown in Fig.[^ 


appears over a pre-defined region in FP2. Then, the resulting 
phase mask needed to provide the necessary phase shift to match 
the optimized PP2 phase is calculated. We note that, for simplic¬ 
ity, all of the simulated optical configurations include the Lyot 
stop, even when there is no VPM or pupil plane mask present. 

For exoplanet imaging, a dark hole in the on-axis PSE in FP2 
is desired starting at a few AjD in angular separation from the 
star location to the edge of the field of view. Here, we chose to 
create a dark hole ranging from about A Aj Dio 20 AjD. The outer 
boundary of the dark hole roughly corresponds to the control 
radius of state-of-the-art extreme adaptive optics systems like 
GPI ( [Macintosh et al.|2014D and SPHERE ( [Fusco et al.|20T^ . 


4, Performance 

During the phase optimization routine, the amount of light that 
appears in the dark region of the on-axis PSE in FP2 decreases 
iteratively. However, the off-axis PSE quality is also affected. 
We find the optimal phase mask by monitoring the contrast and 
throughput at each iteration. The former is defined as the energy 
ratio between on-axis and off-axis PSE, integrated over the full 
width at half maximum (FWHM) of the off-axis PSE. The latter 
is defined as the encircled energy within the FWHM of the off- 
axis PSE, normalized by the encircled energy without the VPM 
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Fig. 5. (a) Contrast and (b) throughput of an off-axis source as a function of angular displacement a for the LPMs optimized to produce a semi- 
annular dark region, (c)-(d) Same as (a)-(b), but for a full annulus. The contrast and throughput afforded by the E-ELT with and without EPMs are 
also shown for comparison. 


and LPM. In other words, the contrast compares the diffracted 
light from an on-axis point source to the signal from an off- 
axis companion of equal magnitude, whereas the throughput is 
a measure of the detected signal from the companion. Here, the 
azimuthal average of the on-axis PSF is used, but only off-axis 
PSFs with displacements along the x-axis are calculated for com¬ 
putational convenience. The algorithm is stopped when the con¬ 
trast improvement between iterations is small, while also limit¬ 
ing the loss in throughput. 

The performance of the VC with an LPM is shown in Fig. 
[^in terms of contrast and throughput of a companion imaged at 
an angular separation a from the star. The values achieved with 
the VC alone, EPMs, and no phase mask (E-ELT) are also shown 
for comparison. The mean contrasts and throughputs within 4-19 
AjD are reported in Table The E-ELT, without a coronagraph 
present, provides a mean contrast of 2.8 x 10“"^. This value is 
reduced to 3.6 x 10“^ with VC2 and 4.9 x 10“^ with VC4. The 
introduction of an LPM signihcantly reduces the contrast to the 
10“^ - 10“^ range for both of the optimization regions consid¬ 
ered. The EPMs offer improvement over the VCs without LPMs, 
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Table 1. Mean contrast and throughput from a - A AjOio a - 19 AID. 


Design 

Contrast 

Throughput 

E-ELT 

2.8 X 10-^ 

1.00 

VC2 

3.6 X 10-5 

0.85 

VC4 

4.9 X 10-5 

0.75 

VC2-hLPMl 

4.0 X 10-^ 

0.46 

VC2-rLPM2 

1.2x 10-® 

0.27 

VC4-rLPM3 

4.2 X 10-’ 

0.43 

VC4-rLPM4 

1.2 X 10“® 

0.26 

EPMl 

2.9 X 10-® 

0.48 

EPM2 

2.1 X 10-5 

0.26 


but yield contrast that is approximately 7-17 times greater than 
the VC-fLPM combination. 

The LPMs and EPMs presented here were designed to yield 
roughly the same throughput. Specihcahy, the throughputs are 
approximately 0.5 for the semi-annular dark region and 0.3 
for the full annulus. Like many advanced coronagraphs, the 
VC-fLPM combinations offer improved contrast performance at 
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the cost of off-axis throughput. In this case, the throughput is 
mainly a function of the inner boundary of the optimization re¬ 
gion, which has been made as small as possible. Specifically, the 
inner boundaries used for VC2 are cr = ?> AjD and cr = 3.25 /l/D, 
for LPMl and LPM3, respectively. For VC4, a = 3.5 AjD 
is used for both LPM2 and LPM4. In the case of the EPMs, 
a = 2.25 AjD and a = 3.5 AjD are used for EPMl and EPM2, 
respectively. We note that significantly better contrast may be 
achieved with the optimization region further from the star. In 
addition, the resulting region of optimum contrast is typically 
smaller than the pre-defined optimization region due to edge ef¬ 
fects. 

In general, the best design of a high-contrast imaging instru¬ 
ment is chosen with several performance aspects in mind, includ¬ 
ing the contrast, size and shape of the discovery region, off-axis 
throughput, and sensitivity to chromatic effects. Each of these 
are addressed in the following discussion. Sensitivity to practi¬ 
cal errors such as imperfect wavefront control, alignment, vibra¬ 
tion, and partial resolution of the star are also important and are 
considered in the next section. 

The desired contrast and dark region shape depends on the 
observational goal; that is, a full annulus is beneficial for disk 
imaging and discovering new companions, whereas only a par¬ 
tial annulus is needed for characterization of a known object. We 
also note that a semi-annulus may be sufficient for 360° high- 
contrast imaging using vector-phase elements and clever polar¬ 
ization tricks (e.g. |Otten et al.|2QT^ . 

For the methods present here, extending the dark region 
within the central lobe of the on-axis PSF does not allow the 
~ 10“^ contrast to be preserved with the VC-i-LPM designs. In 
these cases, the algorithm does not converge to a solution with 
small contrast and large throughput values. Thus, in addition to 
observational goals, the optimal size and shape of the dark region 
depends on the size and shape of the central lobe of the on-axis 
PSF, and ultimately the telescope aperture and the coronagraph 
instrument. We also find that, in comparison to the full annular 
case, the inner boundary of a partial annular dark region may be 
forced closer to the optical axis without significantly degrading 
off-axis throughput performance. 

The outer edge of the optimization region, on the other hand, 
will likely be matched to the control region of the adaptive optics 
system. Here, we optimize within a 20 T/D radius, which is rep¬ 
resentative of current state-of-the-art adaptive optics. We note, 
however, that instruments developed for the future E-ELT may 
have larger control regions thanks to next-generation deformable 
mirrors. Though the outer boundary does not significantly af¬ 
fect the contrast achieved, increasing the outer radius generally 
leads to higher spatial frequency phase variations in the calcu¬ 
lated phase masks. Ultimately, the allowable phase variations are 
limited by the manufacturing process, which will be investigated 
in future work. 

An advantage of the VC-rLPM, over other coronagraph de¬ 
signs such as those with occulting focal plane masks, is its in¬ 
trinsic achromaticity. That is, assuming the phase masks apply 
perfectly achromatic phase shifts, the calculated power leaked 
through the LS, throughput, and contrast do not directly depend 
on wavelength. Rather, the wavelength dependence is limited 
to the scaling of the PSP, which causes radial blurring of the 
dark hole. Polychromatic light results in degraded contrast along 
a narrow annular strip at the edges of the optimization region 
where the affected width is given by roughly half the bandwidth 
fraction (0.5 AT/T) multiplied by the inner and outer boundary 
angle. Por example, with a bandwidth of 20%, the inner and 
outer boundaries of the dark hole are expanded and contracted by 
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Fig. 6. Off-axis PSF’s (in log irradiance) along the v-axis 
for (a) VC2+LPM1, (b) VC2+LPM2, (c) VC4+LPM3, and 
(d) VC4+LPM4 at Of = 10 A/D. 


about 10%, respectively. Achromatic phase masks may be fabri¬ 
cated by direct writing of multilayer li quid crystals ([Komanduri 
|et al.|2013]|Miskiewicz & Escuti|2014||Qtten et al.|2014|). 

The width and shape of the off-axis PSP are very impor¬ 
tant for imaging applications. Pigure shows the off-axis PSP 
for each PPM at a = 10 AjD. A relatively high-quality PSP is 
formed in each case (see Pig.|^b),(d) for corresponding through¬ 
put values) and the shape is maintained for off-axis sources 
throughout the dark region of the on-axis PSP. Moreover, the 
PSP is spatially invariant over the optimization region, which is 
desirable for performing deconvolution on the obtained images 
as well as other forms of post-processing. 

5. Sensitivity to errors 

A primary practical concern is that aberrations may cause 
starlight to appear inside the dark zone and, therefore, degrade 
contrast. In this section, the LPMs are shown to be robust to typ¬ 
ical wavefront error expected for an adaptively-corrected high- 
contrast imaging instrument. The contrast is often limited by 
imperfections in the optical surfaces, which form quasi-static 
speckles in the image. Thus, we model the phase error in the 
pupil with the normalized power spectral density (PSD) function 

where ^ is the magnitude of the spatial frequency and is the 
"cut-off" spatial frequency. In the following simulations, ran¬ 
dom phase screens are generated and scaled to root-mean-square 
wavefront error oj. Por reference, state-of-the-art systems typi¬ 
cally have values = 200 cycles/m and oj ^ 4/100. 

To simulate realistic aberrations, a random phase screen is 
generated at PPl and the resulting contrast is calculated. Pigure 
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co/X 

Fig. 7. Mean contrast for a point source displaced along the x-axis 
(without azimuthal averaging) within a = 4 - 19 A/D sls a function 
of root-mean-square phase error in PPl o) for (a) semi-annular and (b) 
annular dark regions. 


[7] shows the mean contrast for a point source displaced along the 
x-axis within a = 4 - 19 A/D sls a. function of oj. Although un¬ 
wanted speckles appear in the dark region, the VC and LPMs of¬ 
fer contrast improvement in the presence of wavefront errors po¬ 
tentially achieved on current and next-generation high-contrast 
imaging instruments (oj < A/100). Such aberrations levels re¬ 
quire post-coronagraphic low-order wavefront sensing solutions, 
such as those proposed in |Codona & Kenworthy| (|2013|), |Singh| 
etaL] ( [20T^[2to] ), or |Huby et al.| ( |2015| ). 


The sensitivity to low-order aberrations may also be con¬ 
trolled by the choice of focal plane mask. Specifically, the en¬ 
ergy transmitted through the LS increases as for a A/D 
and therefore an / = 4 VPM may be more suitable than / = 2 on 
a system susceptible to pointing errors and/or vibrations. In ad¬ 
dition, reduced sensitivity to tip-tilt avoids leaked light owing to 
partial resolution of the star ( [Delacroix et al.pol^ . We note that 
the EPM designs are the least sensitive to these types of errors. 


Imperfect alignment of the LPM will also negatively effect 
the achieved contrast performance. For the LPMs presented here, 
we expect the mean contrast to increase by a factor of -1.5 for 
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an offset of 0.1% of D. Alignment to within <0.2% is possible 


in practice ( jMontagnier et al.|2007| ). 


6. Conclusions 

We find that phase-only optical elements placed in the Lyot plane 
of a vortex coronagraph may improve the contrast performance 
on telescopes with complicated apertures. The combination of 
a vortex coronagraph and Lyot-plane phase mask provides bet¬ 
ter than 10“^ contrast within a = 4 - 19 A/D on heavily ob¬ 
scured telescopes, such as the E-ELT, at the expense of through¬ 
put. Moreover, we have shown that the improvement offered by 
an LPM is robust to realistic aberrations. Reducing the starlight, 
and its associated noise, enables sensitive high-contrast imag¬ 
ing of circumstellar disks and exoplanets. Phase masks similar 
to those described here are expected to improve the performance 
of current Lyot-style coronagraphs on ground-based telescopes, 
and may provide a route to terrestrial planet imaging with future 
space telescopes. 
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